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Abstract 

The optical excitations of elongated graphene nanoflakes of finite length are investigated 
theoretically through quantum chemistry semi-empirical approaches. The spectra and the re- 
sulting dipole fields are analyzed, accounting in full atomistic details for quantum confinement 
effects, which are crucial in the nanoscale regime. We find that the optical spectra of these 
nanostructures are dominated at low energy by excitations with strong intensity, comprised 
of characteristic coherent combinations of a few single-particle transitions with comparable 
weight. They give rise to stationary collective oscillations of the photoexcited carrier density 
extending throughout the flake, and to a strong dipole and field enhancement. This behavior is 
robust with respect to width and length variations, thus ensuring tunability in a large frequency 
range. The implications for nanoantennas and other nanoplasmonic applications are discussed 
for realistic geometries. 
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In the last few years remarkable interest has grown for nanoplasmonics and the perspectives it 
offers to merge electronics and photonics at the nanoscale.'^ A wide range of potential applica- 
tions can be designed, including sensing and spectroscopic techniques .'SnS fabrication of nanoan- 
tennas'^HH and light emitters ,^2113 ^g ^qW ^g beyond-THz optical devices or solar-energy con- 
version systems.'^^^ While great attention has been devoted to metal nanoparticles, due to the 
relative ease to produce them and to their possibility to support surface modes j^^^HlSl new materials 
and metamaterials are now being exploredjI^SEU which, in addition to enhanced optical responses, 
are able to optimize circuit integration and reduce losses. Graphene has proved to have unique 
electronic and mechanical properties ,221 and has been more recently investigated also for photonics 
and optoelectronics.ESlES] jn the field of (nano)plasmonics, so far most emphasis has been devoted 
to spectroscopy of plasmons in extended graphene, either doped or undopedj ^^ and to the large 
lifetimes of their excitations compared to conventional metals-I^SES Some interesting predictions, 
mostly based on macroscopic models, have been proposed for plasmonics in spatially-modulated 
graphene of micron and sub-micron size range. '22ESEII 

Here we focus on graphene in a completely different regime, typical of nanoscale structures, 
where an optical gap opens as a consequence of quantum confinement, and we show that field 
enhancement effects can also be seen. This regime has become particularly exciting in view of the 
recent production of controlled graphene wires by chemical routes. I^SIHII y^^ consider the case of 
self-standing elongated graphene nanoflakes (GNFs) with H-terminated edges, which we analyze 
by applying a fully-microscopic quantum-chemical approach. These flakes can be thought of as fi- 
nite portions of armchair-edged graphene nanoribbons (GNRs)^-^ [seejTJa)] of variable widths and 
lengths. We focus on the optical properties, which are sensitive to the structural details of the sys- 
tem, namely length and width modulation, by computing their UV-vis spectra. In the low energy 
region intense and tunable peaks are recognized, characterized by a coherent coupling of transi- 
tions that yield "collective" excitations with strong transition dipole. The local field enhancement 
produced by these large dipole excitations is discussed for a prototypical flake of sub-nanometer 
width and nanometer length. 
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Figure 1: (a) Ball-and-stick model of a prototypical elongated graphene flake, characterized by 
length parameter Nl = 28 (x direction) and width parameter Nw = 7 (y direction). Edge atoms 
are passivated with H (white balls), (b) Optical spectrum (dashed line) calculated for the system 
in (a): the oscillator strength (OS) of individual excitations is indicated by black bars. The main 
excitations are labeled according to their polarization with respect to the long axis of the system 
(x), i.e. L for longitudinal and T for transverse polarization, (c) Optical spectra of a series of 
graphene flakes of fixed width (A^iy = 7) and variable length, from 24 to 88 A, labeled after their 
length parameter N^, ranging from 10 to 40. The OS is normalized with respect to the flake length 
in both (b) and (c). All spectra are obtained by introducing a Lorentzian broadening of 25 meV. 



To investigate the nature of these optical excitations, their microscopic origin and their size- 
dependent trends, we adopt the semi-empirical Hartree-Fock-based method ZINDO, which imple- 
ments the configuration-interaction (CI) procedure including single excitations only (S).'^ This 
method is known to provide reliable results for the optical spectra of aromatic moleculesS^BZIHil 
All calculations are performed starting from optimized geometries obtained with AMI ,'^2'^ used 
also for the calculation of mean-field ground state properties. 

The analysis of the optical properties is carried out for a series of prototypical short graphene 
nanoribbons or graphene nanoflakes (GNFs) of fixed width (~ 7.3 A) and variable length (from 
24 to 88 A). Following the standard notation for armchair GNRs,'^ these structures can be char- 
acterized by a width parameter Nw, which indicates the number of dimeric lines along the zigzag 
direction (y axis). In addition, we here introduce a length parameter N^ that corresponds to the 
number of zigzag chains in the armchair direction {x axis), excluding the ends. As shown in the 
model structure of [TJa), the flake edges are passivated with H atoms and the flake ends are shaped 
to minimize the zigzag-edge contributions .1^^^^ The structure of a prototypical flake with Nw = V 
and Nl = 28 is displayed in [TJa) and its UV-vis spectrum (dashed line) is shown in [TJb), where 
the oscillator strengths (OS) of the individual optical excitations are indicated by black bars. The 
optical spectrum is dominated by three intense peaks: the first two, LI and L2, correspond to ex- 
citations with a large transition dipole along the x axis of the system (L stands for "longitudinally 
polarized"); the third peak is found at higher energy and is dominated by the excitation T2, which 
shows a transverse polarization (y direction) with respect to the long axis of the system. In addi- 
tion to these, we also find a transversally polarized excitation Tl with negligible oscillator strength 
below the first peak LI, and a weak longitudinally polarized excitation (LW) between LI and L2. 
[TJ^c) displays the calculated UV-vis spectra for GNFs of the same width (Nw = V) and increasing 
length, with N^ ranging from 10 to 40. The main features described for the case of N^ = 28 are 
maintained along the flake series, except for LI becoming the lowest-energy excitation in place of 
Tl for the longest GNF (Nl = 40): this is a signature of the approaching behavior of quasi- ID arm- 
chair graphene nanoribbons (AGNRs).^-^ Moreover, as a result of confinement, we find an overall 



red-shift of LI and L2, with the optical gap decreasing from 3.14 to 2.48 eV. The energy difference 
between these longitudinal excitations is almost unaffected upon flake length increase, ranging 
from about 0.55 eV (Nl = 40) to 0.60 eV (Nl = 10). As expected, the energy of Tl and T2 does 
not vary considerably with length, since they are polarized along the flake width, kept fixed. At 
increasing length, we also notice a significant transfer of OS to the lowest energy peak LI, which 
tends to saturate for the longest flake, again similarly to the behavior of infinite nanoribbons.'^S 

LI: lowest energy excitation 
(a) transition density 




single-particle states: coherent combination 
(b) (c) 




Figure 2: Lowest energy excitation LI ofN^ = 28 graphene flake shown in[Ija). (a) Transition den- 
sity computed according to ??. (b) Single-particle wave functions mostly contributing to LI. The 
lowest energy excitation LI is obtained from a coherent superposition of single particle transitions, 
as indicated in scheme (c). 



To gain further insight into the nature of these spectral features, we analyze the composition 
of all low-energy excitations in terms of molecular orbital (MO) transitions and CI weights. The 
numerical details are presented in the Supporting Information (Table SI) for few selected cases; 
here we focus again on Nl = 28, starting from the first bright excitation LI. As illustrated in|2} 
this excitation arises from the coherent combination of single -particle transitions involving several 



harmonics of the same states. Note that, while in shorter flakes the HOMO — )■ LUMO transition 
contributes for most of the weight, in longer structures several higher harmonics enter the compo- 
sition with comparable weights. The fact that transitions between occupied and virtual MOs with 
the same envelope function modulation concur to form the excitation helps explaining the large OS 
observed, as this mechanism tends to maximize the wave function overlap. The spatial extension 
of the MOs increases with their energy distance from the frontier orbitals, further contributing to 
the large transition dipole. This effect is made particularly evident by plotting the transition density 
of the excitation, defined as: 

ia 

where Cf^ are the CI coefficients of the ;?-th excited state, corresponding to single excitations from 
the occupied ^,- to virtual (pa orbitals.!^ From [2[a) we notice that the density extends homoge- 
neously over the entire structure with an underlying dipolar character, and is not concentrated just 
to the central portion. This is evident also by inspecting the contour plot of the corresponding 
electric potential, pictorially represented in|3j 
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Figure 3: Pictorial view of the electric potential generated by the transition density relative to the 
excitations LI, L2, LW and T2 of the Nl = 28 flake. 



Consistently with their polarization, LI and L2 have a large dipole oriented along x, while that 
of T2 is along y. As shown for LI in [2} also L2 and T2 arise from combinations of transitions 
between states with the same envelope function modulation. On the contrary, LW is mainly com- 
posed by transitions involving MOs with different number of nodes. This gives rise to a multipolar 
modulation in the electric potential, as shown in|3} and leads to the weaker intensity observed for 
this excitation. 
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Figure 4: Field enhancement of longitudinally polarized excitations LI, L2 and LW of Nw=^ and 
Nl=2S graphene nanoflake computed along the flake axis, at 3.3 A from its plane. Close to the flake 
center, above its surface, the field enhancement is susceptible of the transition density modulation 
[see[2]^a)]. In the inset, pictorial sketch of the field enhancement contour plot in {x,y) and {x,z) 
views. Notice that the values of the field enhancement are sensitive to the choice of Tp, according 
to ??: here a conservative value rp=25 meV is adopted. 

The large dipole strength and the "beyond-single-particle" nature of LI and L2 recall the main 
features related to collective plasmonic excitations in small metal clusters, as analyzed by means 
of atomistic methods. I^SHSH-wTe ^j^^^g calculate for our systems a quantity which is usually adopted 
for the characterization of plasmonic excitations, i.e. the local enhancement of an electromagnetic 
field incident on the system. To this end, we approximate the flake response function, at resonance 
conditions, via the electronic density variation dp{(OQp\r) induced by the periodic external electric 
field EQexp{—iG)Qpt), given by: 






(2) 



Here ^l^P is the transition dipole for the excitation to the in-resonance p state, and Tp is the cor- 
responding decay rate, related to the intrinsic absorption bandwidth. Details on the derivation of 
?? are given in the Supporting Information. The oscillating charge distribution dp{o)op',r) orig- 
inates an additional electric field E e(r)exp{—icoopt) leading to an overall enhanced electric field 
(Eq + Ee) nearby the graphene flake. The transition dipole jU*^^ is a straight result from the cal- 
culation, while Tp is a critical parameter in determining the maximum enhancement [the larger 
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this value, the smaller the enhancement, as shown in ??]. The field enhancement is obtained by 
deriving the electric Coulomb potential generated by the transition density, and expressed in units 
of £0 (Ee/Eo). Since experimental data are not yet available for graphene nanoflakes such as those 
examined here, we adopt in the figures a conservative value rp=25 meV, chosen from the spectral 
hnewidths measured for single-wall carbon nanotubes (SWCNTs) at room temperature in solu- 
tion.SSlDue to the critical dependence of the field enhancement on the choice of Tp [see ??], it is 
worth noting that the field enhancement shown in |4] would increase of over one order of magnitude 
for Fp ~ 1 meV, as found for suspended SWCNTs at low temperature.'^ In H we show the field 
enhancement produced by the three main longitudinal excitations (LI, LW and L2) identified in 
[ijb) for the GNF of Nw=^ and Nl=2S, computed along the flake longitudinal (x) axis, at a distance 
of 3.3 A from the basal plane [see[5[b)]. The oscillating character of the curves in the region above 
the flake surface is due to the transition density modulation [see|2];a)]. Beyond the flake border 
the oscillations disappear and, at distance much larger than the flake half-length (~30 A), the field 
enhancement assumes a Coulomb-like decay, proportional to 1/r^. According to the OS of the 
corresponding excitations, the field enhancement for LI and L2 have comparable values, while 
that of LW is about three times lower. 

Finally, we discuss the stability and the tunability of these optical properties with respect to 
length and width modulation. Focusing specifically on the lowest energy peak LI, the computed 
values of field enhancement are basically independent of the ribbon length (see Supporting In- 
formation, Figure SI). On the other hand it is well known that both the electronic and optical 
properties of quasi- ID AGNRs are sensitive to width modulation.l ^^ ' ^^" ^ Also for the finite flakes 
considered here three families are identified,'^ characterized by different electronic gaps, with the 
smallest values pertaining to the Nw = 3p + 2 family (p integer). Here we investigate two addi- 
tional graphene flakes of width parameters Nw=6 and N\y=^, keeping their length fixed at Ni=2S. 
As shown in [5]; a), the three main peaks (LI, L2 and T2) observed for Nw=^ are preserved in all 
families. While T2 redshifts at increasing width, as expected, the energy and the intensity of the 
longitudinal excitations LI and L2 are closely related to the electronic properties of each family. 
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Figure 5: (a) UV-vis spectra of a series of graphene flake of fixed length (Ni = 28) and variable 
width, with Nw ranging from 6 to 8. A Lorentzian broadening of 25 meV is introduced, (b) Field 
enhancement of LI excitations along the flake longitudinal axis, at a distance of 3.3 A from its 
basal plane, according to the scheme reported on top of the graph. Notice that the values of the 
field enhancement are sensitive to the choice of Tp, according to ??: here a conservative value 
r„=25 meV is adopted. 



10 



In particular the trend for the electronic gaps^^ is reflected in the excitation energies of LI. The 
field enhancement computed for the LI excitation of each GNF of different width is shown in[5]^b): 
we find a distinctive enhancement for A^iy=8 compared to the other two, which already comes from 
the transition dipole moments of the individual single-particle excitations. 

In summary, we have analyzed the optical excitations of finite graphene nanoribbons, of sub- 
nanometer width and nanometer length. At low energy the UV-vis spectra are dominated by intense 
excitations with longitudinal polarization with respect to the flake long axis. These are character- 
ized by a "collective" character, coming from coherent superposition of MO transitions with the 
same envelope function modulation. The investigated excitations are tunable both in energy and 
intensity upon appropriate length and width ribbon modulation. The field enhancement computed 
for these excitations suggests the applicability of these systems as nanoantennas and in other opto- 
electronic and nanoplasmonic applications. 
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Excitation Analysis in Terms of Orbital Transitions and Weights 

In Table SI we report the lowest energy excitations, labeled after Figure 1(b) in the main text, of 
three graphene nanoflakes (GNFs) at increasing length. We consider the shortest (A^l=10, ~ 24 A) 
and the longest (A^^=40, ~ 88 A) ones, in addition to Nl=2S GNF (length ~ 63 A) discussed in 
details in the main text. The excitation energy, oscillator strength (OS) and composition in terms 
of molecular orbital (MO) transitions for the main excitations are presented. 

From the results presented in Table SI, the gain of OS characterizing the first bright excitation 
LI at increasing length is evident. It is also worth noting that for the longest considered ribbon 
(A^L=40) the lowest energy excitation is indeed LI, being basically degenerate with Tl. 



System 


Excitation 


Energy [eV] 


OS 


Transitions (weight) 




Tl 


2.54 


0.0005 


HOMO-1 -^ LUMO (0.37) 
HOMO -^ LUMO+1 (0.39) 




LI 


2.72 


1.62 


HOMO -^ LUMO (0.79) 


A^L=10 


L2 


3.32 


5.43 


HOMO-3 -^ LUMO+3 (0.12) 
HOMO-1 -^ LUMO+1 (0.76) 


LW 


3.48 


0.34 


HOMO-4 -^ LUMO (0.14) 

HOMO-2 -^ LUMO+2 (0.26) 

HOMO -^ LUMO+4 (0.31) 




T2 


3.59 


1.56 


HOMO-1 -^ LUMO (0.41) 
HOMO -^ LUMO-i-1 (0.39) 



Continued on next page 



S2 



Table SI - continued from previous page 



System 


Excitation 


Energy [eV] 


OS 


Transitions (weight) 




Tl 


2.48 


0.00 


HOMO-3 ^ LUMO (0.21) 
HOMO ^ LUMO+3 (0.22) 










HOMO-2 -^ LUMO+2 (0.11) 




LI 


2.51 


8.19 


HOMO-1 -^ LUMO+1 (0.21) 
HOMO -^ LUMO (0.48) 










HOMO-2 ^ LUMO (0.21) 




LW 


2.78 


0.35 


HOMO-1 ^ LUMO-i-1 (0.15) 


Nl=IA 








HOMO ^ LUMO-i-2 (0.23) 










HOMO-6 -^ LUMO-i-6 (0.12) 




L2 


3.07 


10.00 


HOMO-4 -^ LUMO+4 (0.20) 
HOMO-3 -^ LUMO-i-3 (0.48) 










HOMO-4 ^ LUMO-i-2 (0.16) 




T2 


3.57 


3.35 


HOMO-3^ LUMO (0.18) 

HOMO-1 -^ LUMO-i-4 (0.15) 

HOMO -> LUMO-i-3 (0.17) 
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Table SI - continued from previous page 



System 


Excitation 


Energy [eV] 


OS 


Transitions (weight) 










HOMO-2 -^ LUMO+2 (0.13) 




LI 


2.48 


12.56 


HOMO-1 -^ LUMO+1 (0.19) 
HOMO -^ LUMO (0.36) 




Tl 


2.48 


0.00 


HOMO-4 -^ LUMO (0.15) 
HOMO -^ LUMO+4 (0.16) 










HOMO-2^ LUMO (0.18) 


A^L=40 


LW 


2.68 


0.62 


HOMO-1 -^ LUMO+1 (0.11) 
HOMO -^ LUMO+3 (0.19) 










HOMO-7 -^ LUMO+7 (0.12) 




L2 


3.03 


13.30 


HOMO-5 -^ LUMO+5 (0.19) 
HOMO-4 -^ LUMO-i-4 (0.35) 




T2 


3.59 


5.16 


HOMO-4 -^ LUMO (0.14) 
HOMO -^ LUMO-i-4 (0.13) 



Table SI: Energy, oscillator strength (OS) and composition of the main excitations characterizing 
the spectra of A^iy=7 GNFs of length A^l=10, Nl=2S and Nl=40. In the last column we include the 
MO transitions with relative weights larger than 0.1. 

In Table S2 we report the lowest energy excitations of A^vK=6, A^iy=7 and Nw=^ graphene flakes 
of length Ni=l'&. Excitation LI in the largest flake (A^vK=8) presents an increased OS of about 10% 
with respect to Nw=6 and Nw=^ GNF, where it keeps almost the same intensity. Also the energy 
of LI is red shifted of over 1 eV in A^iy=8 GNF, with respect to that of the other two, as observed 
in the UV-vis spectra of Figure 5(a) of the main text. 
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Table S2: Energy, oscillator strength (OS) and composition of the two lowest energy excitations of 
graphene nanoflakes with Nl=2S and different width, namely A^w'=6, Nw=^ and Nw=^- In the last 
column only we include the MO transitions with relative weights larger than 0.1. 



System 


Excitation 


Energy [eV] 


OS 


Transitions (weight) 


Nw=6 


LI 


2.47 


8.23 


HOMO-1 -^ LUMO+1 (0.20) 
HOMO -^ LUMO (0.54) 




Tl 


2.48 


0.002 


HOMO-3 -^ LUMO (0.21) 
HOMO -^ LUMO+3 (0.22) 


Nw=l 


LI 


2.51 


8.19 


HOMO-2 ^ LUMO+2 (0.11) 

HOMO-1 ^LUMO+1 (0.21) 

HOMO ^ LUMO (0.48) 


A^w=8 


LI 


1.37 


9.15 


HOMO-1 ^ LUMO-i-1 (0.18) 
HOMO -^ LUMO (0.69) 



(a) 



(b) 
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Figure SI: Field enhancement of longitudinally polarized excitation LI for selected graphene 
nanoflakes of fixed width Nyv=l and increasing length from about 24 {Nl=10) to about 75 A 
{Nl=3A), computed along the flake axis (a) at 3.3 A from its plane and (b) along the ribbon plane, 
starting from its lateral end. Notice that the values of the field enhancement are sensitive to the 
choice of Fp, according to Eq. (2) in the main text: here a conservative value rp=25 meV is 
adopted. 
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Effects of Length Increase on the Field Enhancement 

We report in Figure SI the field enhancement produced by LI excitation along the longitudinal 
(x) axis of selected graphene flakes of increasing length, at different heights from the basal plane: 
at 3.3 A above the plane from the flake center [see Figure Sl(a)] and along the plane from the 
flake end [see Figure Sl(b)]. It is worth noting that above the ribbon plane the field enhancement 
is basically independent of the length, and the curves are characterized by an oscillating behavior 
related to the transition density modulation. On the other hand, further from the flake lateral 
border, where the dipolar character of the transition density becomes predominant [see also Figure 
2(a) in the main text], the field enhancement curves follow the oscillator strength gain of LI upon 
increasing length [see Table SI and also Figure 1(c) in the main text]. 

Enhanced Response Field from Transition Density 

We derive the relation between the enhanced field and the transition properties of the molecule. 
In general, the electric field Eg produced by the polarization of the molecule due to an external 
potential Vq oscillating at a frequency ft) (for example, due to an external electric field Eq considered 
in the quasi-static limit) can be found from: 

V-Ee = 4-7J:dp{f) (1) 

where 5p{r) is the electronic density induced in the position r by Eq. From response theory, it 
comes: ^ 

5p{a);r) = Jxio};lr)Vo{r)dr', (2) 

where the molecular response function can be written as: 

y(^.p/N y Poiir)pd7') PiQ{7)pm{r') 

^^ ' ' ^ ^Qhco + ihTi + Ei-Eo -hco-ihTi + Ei-Eo ^ 
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In Eq. (3), po/ are the transition densities between the ground state and the excited state i, and 
Ei are the energies of the states. When the frequency ft) is close to a molecular transition from 
to p (i.e., ft) ^ (Opo = {Ep — Eo)/h), one denominator in Eq. (3) becomes very small, and the 
corresponding term dominates the sum over /, so that: 

z(«„o;^/)^ P°'-'2^;'"" . (4) 

Therefore: 

Sp{co;r) = P^ J pp^(r')Voir) dr' . (5) 

If we now assume that VQ{r') is the electrostatic potential associated with an external probing field 
in the quasi static limit, then: 

VQ{f') = -f'-EQ (6) 

and 

5p(ft>op;r) = -^^AlpO-^o, (7) 

where jUpo is the molecular transition dipole, evaluated here by the ZINDO/S method. By solving 
Eq. (1) using the expression of 5p(ft);r) in Eq. (7), the electric potential and the field Eg plotted in 
Figure 3 and Figure 4 of the main text are finally obtained. 
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